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Abstract. We extend Sklyanin's method of separation of variables to quantum inte- 
grable models associated to elliptic curves. After reviewing the differential case, the 
0> ' elliptic Gaudin model studied by Enriquez, Feigin and Rubtsov, we consider the differ- 

ence case and find a class of transfer matrices whose eigenvalue problem can be solved 
by separation of variables. These transfer matrices are associated to representations of 
, the elliptic quantum group Er^nish) by difference operators. One model of statistical 

mechanics to which this method applies is the IRF model with antiperiodic boundary 
, conditions. The eigenvalues of the transfer matrix are given as solutions of a system of 

' quadratic equations in a space of higher order theta functions. 



^ ! 1. Introduction 

^ I The method of separation of variables in integrable lattice models, proposed by 

O ' Sklyanin, is a method to find eigenvalues and eigenvectors of transfer matrices. It is 
^ . an alternative to the Bethe ansatz and works in some situations where the Bethe ansatz 
Q>^ ! does not, and gives an insight in the completeness of the Bethe ansatz. The method is 
^ I closely related to Baxter's method (Chapter 9 of 0), and in fact the eigenvalue prob- 
^ ■ lem in the separated variables (in the difference case) becomes the Baxter difference 
. equation. In the Gaudin model, one of the simplest quantum integrable systems, this 
method relates the problem of finding common eigenvectors of Gaudin Hamiltonians to 
the problem of finding differential equations on the Riemann sphere with regular singu- 
lar points whose monodromy is trivial. As noticed by Feigin and Frenkel (see 0), this is 
a special case of the Beilinson-Drinfeld "geometric Langlands correspondance" relating 
I certain local systems on a complex curve to P-modules on moduli spaces of principal 
bundles on the curve. 

Both for quantum integrable models and for the connection to the Langlands program, 
it may be interesting to extend the method of separation of variables to more general 
models. 

The class of quantum integrable systems (families of commuting operators) one con- 
siders in this context arise in different classes. There are the differential models, such 
as the quantization of the Hitchin systems. They are given by families of commuting 
differential operators and are associated to complex curves. The Gaudin operators are 
the operators associated to the Riemann sphere. More generally one considers difference 
or g-deformed models, such as the six-vertex model, which degenerate to the differential 
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model when the parameter q tends to one. These models appear in three sorts: ratio- 
nal, trigonometric and elliptic, depending on the type of coefficients in the commuting 
operators. 

From the point of view of representation theory, the differential models are related to 
Kac-Moody Lie algebras, and the g-deformed models to (infinite dimensional) quantum 
groups. 

We consider here models related to s/2. In the differential rational (genus zero) case 
the separation of variables was considered by Sklyanin |]T0| and Frenkel 0. A version of 
the separation of variables for the genus one differential case was considered by Enriquez, 
Feigin and Rubtsov |Q, who also made an explicit connection to the Langlands corre- 
spondance. The g-deformed rational and trigonometric case were studied by Sklyanin 
|T0| and Tarasov-Varchenko |12|. The latter authors introduce the notion of a difference 



equation with regular singular points, thus giving a g-version of the relation described 
above for Gaudin models. 

Here we consider the g-deformed elliptic case. The class of difference operator we give 
is both a g-deformation of the Enriquez-Feigin-Rubtsov differential operators and an 
elliptic version of the operators studied by Tarasov-Varchenko. Common eigenfuncions 
may be constructed by the method of separation of variables. Moreover the commuting 
difference operators we introduce can be restricted to functions on a finite subset of 
points. This restriction turns out to give the commuting transfer matrices of interaction- 
round-a-face models with antiperiodic boundary condition. This model provides an 
example of a model solvable by separation of variables but not by Bethe ansatz. Other 
elliptic models, related to the XYZ model, have recently been studied by the method of 
separation of variables by Sklyanin and Takebe [|ll] , 



The algebraic structure at the origin of our constructions is the elliptic quantum 
group Er,>ri{sl2). Indeed, the starting point is the construction of a representation of 
this quantum group by difference operators which generalizes the "universal evaluation 
module" of 0. 

The paper is organized as follows: 

In Section 2 we review the separation of variables in the differential elliptic case. Most 
of this part is essentially taken from but we add some remarks on the Bethe ansatz 
and its completeness. 

In Section 3 we explain what is needed from the theory of elliptic quantum groups and 
introduce a class of representations of Erri{sl2) by difference operators and relate them to 
known representations. Twisted commuting transfer matrices are then introduced and 
the method of separation of variables is applied to construct (Bethe ansatz) eigenvectors. 

In Section 4 we consider the restriction of the transfer matrix associated to the tensor 
product of n fundamental representations to functions of a finite set of cardinality 2", 
and show that we obtain the transfer matrix of an IRF model with antiperiodic boundary 
conditions. The eigenvalues are then obtained as the solutions of a system of n quadratic 
equations in an n-dimensional space of theta functions of order n. 
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In an appendix we give an account on "elliptic polynomials" , which are (twisted) theta 
functions of order n. 



2. The differential case 

Let us start by introducing a family of commuting differential operators associated to 
an elliptic curve with n marked points and n highest weight representations of s/2(C). 

Let the elliptic curve he E = <C/T with F = Z + rZ and Imr > 0. The marked 
points are the projections of n points Zi, . . . ^z^ G C, with Zi ^ Zj mod F for i 7^ j. 
The representations are Ma^ , • • • , Ma„ where, for A G C, Ma denotes the Verma module 
(defined below in ^?T]of highest weight A of s/2(C). 

Thus our parameters are r, 21, . . . , Ai, . . . , A„. 

Let d{z) = - EjGzexp(z7r(j + 1/2)V + 27ii{j + l/2){z + 1/2)) be the odd Jacobi 
theta function, and set (T\{z) = • unique meromorphic function of z 

regular on C — F, with a simple pole with residue one at 0, and such that a\{z + r + sT) = 
e^^'^'^axiz), r,s eZ. 

Let e = (00), / = (5o), h = (oi^) be the standard generators of s/2(C). For a G sZ2(C) 
let a^*) denote the action of a on the ith factor of the tensor product M = Maj®- ■ ■(S'Ma^. 
Introduce the following endomorphisms of M depending on 2;, A G C: 

1=1 ^ ^' i=l 1=1 

The family of commuting differential operators is then obtained by the following gener- 
ating function, which is an elliptic version of the generating function of Gaudin Hamil- 



tonians pO[] 



Theorem 2.1. Let for z G C, S{z) he the differential operator acting on functions of one 
complex variables A with values in the zero weight space M[0] = {f G M, h^'^^v = 0} 
ofM. 



d I " ' 



S{z) =[jy^- ^Hz)j + e,{z)f,{z) + f,{z)e,{z). 

Then S{z)S{w) = S{w)S{z). 

One way of proving this theorem is to notice that it is a special case of the flatness 
of the KZB connection (Prop. 2 in [Q). The relation to the KZB connection is the 
following. In the sZ2(C) case, the KZB connection involves the differential operators 
(appearing on the right hand side of the KZB equations) 



.t^. 2 0{zj - Zk 
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j = 1, . . . ,n and 



(9A2 2 ^ 2 ^(2. - 2fc) ' ' ' ' (9A 

j,k=l 

The fact that the connection is flat means in particular that these operators form a 
commuting family when acting on M[0]- valued functions. The terms with j = k are 
understood as the limit as the argument zj — tends to zero. So 9"{{))/9{{)) means 
e'"(0)/^'(0) and {dax/d\){Q) = {9'/9y{\). Note that EJ=i vanishes on M[0]-valued 
functions. We call the operators Hj, j > 1, the elliptic Gaudin Hamiltonian and Hq the 
(generalized) Lame Hamiltonian (it is the Lame operator if n = 1). 
The relation between these commuting operators and S{z) is: 

Proposition 2.2. Le^ C{z) = 9'{z)/9{z), p{z) = -('(z). Then 

n n 

'^(^) = 5Z y^*^^ ~ + ^HkC{,z - Zk) + Ho, 

k=l k=l 

and Cj is the Casimir value Cj = |Aj(Aj + 2). 

The proof of this fact follows by noting that S{z) is a meromorphic doubly periodic 
function of z with at most double poles at the points zj. By expanding S{z) in a 
Laurent series up to the constant term at the z = Zj, we find that the difference between 
left-hand side and right-hand side is a differential operator whose coefficients are regular 
elliptic functions vanishing at least at one point. Such an operator vanishes by Liouville's 
theorem. 

The eigenvalue problem for common eigenfunctions of Hq, . . . , Hn can then be formu- 
lated as S{z)u = q{z)u with q{z) = Yl'k=i ^Pi^ - ^k) + YIk=i ^kCi.^ - Zk) + eo. The 
eigenvalue of Hj is then e^, j = 0, . . . , n. Since Ylij>i — 0' must necessarily have 
ei + ■ ■ ■ + e„ = 0. 

Common eigenfunctions of Hj and thus of 5* (2;) can be obtained by the Bethe ansatz 
method. They have the form f{wi) ■ ■ ■ f{wm)vo where vq is the tensor product of highest 
weight vectors, ""^ = | X] Wi, . . . , Wm are a solution to the system of Bethe ansatz 

equations, see ^ and below. 

2.1. Separation of variables. We realize the representations of s/2(C) by differential 
operators: 

Lemma 2.3. For any A G C, the map f ^ t, + A, e t-^ + ^Tt defines 

a representation of sl2{C) on C[t], the Verma module M\. If A is a nonnegative integer, 
t^+^C[t] is an invariant subspace and the quotient La = C[t]/t^^^C[t] is irreducible with 
highest weight vector 1 G Ker(e) of weight A. 

^The relation of these functions with the classical Weierstrass ( and p functions (^{z) = ^ + 
where 2771 = 0"'(O)/36''(O). 
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The proof consists of checking the relations [e, /] = h, [h, e] = 2e, [h, f] — —2/ and 
et^+i = 0. 

Therefore we may reahze the tensor product M as C[ti, . . . ,tn\, and the tensor prod- 
uct of irreducible representations (for integer A^) as C[ti, . . . , ^^.(t^^'''^C[ii, . . . , 
Then M[0] consists of homogeneous polynomials of degree m = ^ Afc/2. We may then 
view e\{z), f\{z), h{z), S{z), Hj as differential operators in n + 1 variables X,ti, . . . ,tn- 
They commute when acting on functions which are homogeneous in ti, . . . , t„ of degree 
m. 

2.2. Separated variables. We express the differential operators S{z) in terms of new 
variables so that the eigenvalue problem is reduced to ordinary differential equations. 
Following Sklyanin's idea, the new variables C, yi, . . . , y„ are the zeros and the leading 
coefficient of the operator fx' 

Since fx{z) is realized as a multiplication operator and both sides of this equation are 
functions of z with definite transformation properties under translations by the lattice 
r, this equation does define, locally around a generic point, a biholomorphic change of 
variables (C, yi, . . . , yn) ^ (A, ti, . . . , The formulae are 



U = C- 



n, 0{z. - yj) 



and 

n 

Prom these formulae we deduce the transformation properties of partial derivatives: 

_d_^d_ ^ e'{yj - Zk) ^ d 

The next step is to remark that a function u{C,yi, . . . obeys S{z)u — q{z)u with 
= ELi ^Pi- - ^k) + ELi^fcC(^ - Zk) + eo and Y.j>i^3 = if and only if it 
obeys S{yj)u = q{yj)u for all j = 1, . . . , n and all generic points (i/i, . . . , yn)- Here the 
ambiguous notation S{iij) means: write differential operator S{z) with the coefficients 
on the left of the partial derivatives and replace z by yj in the coefficients. To prove this 
statement notice that, if S{yj)v — q{yj)v, then n^(-^ ~ ^di^i^) ~ ■ ■ ■ ^Z/n) is 



a holomorphic theta function in z of order n vanishing at n generic points %){. It thus 
vanishes, see the Appendix. 

It is then convenient to use the identity \e\{z), f\{z)\ = —h\z) on M[0], to write S{z) 
as S{z) = — ^h{z)Y — h'{z) + 2fx{z)ex{z). so that the last term vanishes if we set 
z = Uj and we get 

Proposition 2.4. A function u{X,ti, . . . ,tn), homogeneous of degree m = \ Yli^k in 
the ti, is a local solution of the partial differential equations S{z)u = q{z)u, z & C if and 
only if 

u{X,ti, . . . ,t„) = Cviyi, ...,yn) 

and V obeys 

(2) = qiyj)v, V,,, = ^ " E f^^^^^' ~ 

k=i 

2.3. Interpolation formula. The formula (|l]) expresses the values at . . . ,?/„ of the 
coefficients of S{z) for z = yj. Since the coefficients of S{z) — | Ckp{z — Zk) are elliptic 
functions of z with at most simple poles at zi, . . . ,Zn, they are uniquely determined by 
these values, and can be calculated by an interpolation formula: let us write S{z) in the 
form 

^ n n 

S{z) = - 'Y^Ckp{z - Zk) + WO{,z- Zk)~^S{z), 

k=l k=l 

SO that S{z) is a theta function of order n. Thus (see the Appendix) 



11 



with 



k=l \\ ^ k=l / k=l 

2.4. Bethe ansatz. The separated equation reads 

n n 
k=l k=l 



Zk 
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with = ^ — X]fc=i ^({y ~ ^k)- It is a second order ordinary differential equation 
with regular singular points at Zk and characteristic exponents and + 1. Following 
Hermite's method to solve the Lame equation (see |T^) we seek solutions of the form 



(4) viy) = e'^yl[eiy-Wk). 

k=l 

Functions of this form are called elliptic polynomials, see the Appendix. 

Let us first assume that Wk 7^ zj mod F for all j, k. Then we also have Wk 7^ Wi for 
k ^ I, since the only solution vanishing with its derivative at a regular point is the trivial 
solution. Rewrite the equation in the form v"{y) — J2k-^kC{y ~ ^k)v'{y) + h{y)v{y) = 
so that b{y) has at most simple poles at the zj. Taking derivatives of a function v of the 
form (^) and setting y equal to one of its zeros, we find the relation: 

v"iwk) = ^C(^^fc - Wj)v'{wk) + 2cv'{wk). 

Inserting this into the differential equation, we see that f is a solution if and only if its 
zeros Wj obey the "Bethe ansatz equations" 



^ ^lC,{wj -zi)- ^ C{wj - Wk) = 2c, j = 1, 

/=1 k:kf!=j 



m. 



Let us now consider the more general case of elliptic polynomials (1) vanishing at zi, for 
i in some subset / of {1, . . . , n}. Since the characteristic exponents at Zi are 0, Aj + 1, a 
solution vanishing at Zi must vanish to order Aj + 1 and is thus divisible by 0{y — 2^)^*^^ 
(this is only possible if Aj G Z>o). Then v{y) = Yliei^iv ~ Zi)~^^^^v{ii) is again of the 
form (H), but with m replaced by m = m — X]je/(^* + It obeys the equation (^ with 
Aj replaced by — Aj — 2 for i e /. 

Thus all elliptic polynomial solutions of (0) are of the form 

rh 

v{y) = n Oiy - z.;)'"^^^' n - ^^)' 

ie/ k=l 

for some subset / of {j \ Kj G Z>o}, such that Wj ^ Wi ^ Zi (j 7^ /) and Wi, . . . , w^, c 
obey the Bethe ansatz equations 



(5) ^ hC,{wj - zi)- ^ C{wj - Wk) = 2c, 

1=1 k:k^j 



m. 



Here A; = —A; — 2 if / G / and A^ = A; otherwise. 

To each such solution there corresponds a common eigenfunction u which, expressed 
in the separated variables, is m = Y\ v{yi). Up to a nonzero constant we get 

(6) u{X) = e''fM---f{wrn)vi. 
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Here Vj = Y[ieiif^^^)^^~^^^o is a product of singular vectors. Only eigenvectors corre- 
sponding to / = have non-trivial projections to eigenvectors with values in the tensor 
product of irreducible representations. 

Eigenvectors of the form (P), such that wi, . . . ,Wm,c are a solution of the Bethe 
ansatz equations (||) with Wj ^ Wk niod F, (j 7^ k) and Wj 7^ Zi mod V are called 
Bethe eigenvectors. 

2.5. Completeness of Bethe eigenvectors. Let us consider the common eigenvalue 
problem 

(7) Hiu{X) = eiu{X), i = 0,...,n, 

A natural class of functions preserved by the operators Hi is given by meromorphic 
sections of a flat line bundle on E. Namely, let for a character x : F ^ C^, 7Y(x) be the 
space of meromorphic functions A 1-^ u{X) G M[0] such that u{X + 1) = x(1)m(A) and 

u{X + t) = xiT)e'''^^'^^^'\{X). 

It is easy to see that ex{z), f\{z) and d/dX — h{z)/2 preserve functions with these 
transformation properties, so that S{z) and Hj preserve 7i(x)- 

It is then natural to look for eigenfunctions (non-trivial solutions of the differential 
equations (0)) in TC{x)- 

Let be the set of (eo, . . . , e„) G C"+^ such that there exists a non-trivial function 
u G 7Y(x) with HjU = eju, j = 0, . . . , n. Let a G F* be the character such that 
a{r + st) = (-1)^+^ 

Theorem 2.5. Let x ^ T*- Then (eo, . . . , e„) G S(x) if and only if J2k>i^i: ~ ^ ^'^^ 
the separated problem 

^^(2/) ~\Y1 ^'^(^'^ + - ^kHy) = (J2 ((y - ^k)ek + eo)v{y), 

k k 

admits a non-trivial elliptic polynomial solution v G 0-m(c"'"x)- -^'^ this case there is a 
Bethe eigenvector with these eigenvalues. 

Proof: A common eigenfunction, viewed as a polynomial in ti has the form 

u{X,t^,...,tn)= ^-i,...,m„(A) JJC'. 

TiiH \-m„=m i 

Replacing tj, A as functions of the new variables, we obtain 

u{X,ti, . . . = Cviyi, . . .,yn), 

where v{yi, . . . , ?/„) is a linear combination of products of theta functions in each of the 
yj with coefficients Umi,...,mn{Yl Vj ~ J2 ^k) and v{yi, . . . , yn) obeys, in each variable, the 
separated second order equation (||). A priori the meromorphic function v{yi, . . . ,yn) 
may have poles on the hyperplane YlVi ~ Yl^i Hiod F. However this is impossible: 
consider f as a function of, say yi with the other variable fixed at some generic position. 
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Then v, as a function of yi, being a solution of a linear second order equation may only 
have singularities at the poles zj of the coefficients. Moreover, since u is in H{x), the 
functions yi v{yi, . . . ,yn) belong to 6m(cr™x)- Thus the separated problem admits a 
non-trivial solution in 0m(c"™x)- As shown in the previous subsection, such a solution 
gives rise to a Bethe eigenvector. □ 

3. The difference case 

3.1. Representations of the elliptic quantum group £'t-^(s/2)- The difference ver- 
sion of the differential operators ex{z), fx{z),dx — h{z)/2, are operators obeying the 
relations of the elliptic quantum group £'^^^(5/2)- Let us recall the definitions 0: we fix 
two complex parameters r, 77, such that Im(r) > 0. The definition of £^,-^^(5/2) is based 
on a dynamical i?-matrix R{z, A) which we now introduce. Let 

^^(A + 2rM^ eiX + z)ei2r^) 

- e{mz-2r^y - ~ e{X)9{z - 2r^y 

Let be a two dimensional complex vector space with basis e[l], e[— 1], and let Eije[k] = 
6jke[i], h = Ell — E_i _i. Then, for 2;, A G C, R{z, A) G End(K ® V) is the matrix 

R{z,\) = Eii0 E11 + E_i^_i0 E_i^_i + a{z,\)Eii0 E^i^^i 

+ a{z, -A)£_i,_i ® £11 + p{z, A)Ei,_i ® E^i^i + p{z, -A)£_i,i ® £i,_i. 

This i?-matrix obeys the dynamical quantum Yang-Baxter equation 

R^^^\z-w, A- 277/1(3) X) i?(23)(^^ ^_2r//i«) = 

{z-w, A) 

in End{V i^V 0V) , z,w, \ E C. The meaning of this notation is the following: i?(^^)(A — 
2T]h^^^) f 1 (g> f 2 ® "^3 is defined as 

{R{z, A - 277/23) f 1 (g) ^2) ® t;3, 

if /if 3 = /i3f3. The other terms are defined similarly: in general, let Vi, . . . , be modules 
over the one dimensional Lie algebra 1) = Ch with one generator h, such that, for all i, Vi 
is the direct sum of finite dimensional eigenspaces Vi[/i] of h, labeled by the eigenvalue 
fi. We call such modules diagonalizable f) -modules. If X G End(l^) we denote by 
G End(ViC>?>- ■ -^Vn) the operator ■ ■ ■(g'ldC^X^Id®- • ■ acting non-trivially on the ith 
factor, and if X = ^ X^on G End{Vi(^Vj) we set X^*-') = ^ X^'^Y^^l If X(/ii, . . . , 
is a function with values in End(Vi ■ ■ ■ ® Ki), then X{h^^\ . . . , h^"'^)v = X(/ii, . . . , fin)v 
if /i^*)?; = /ijW, for alH = 1, . . . , n. 

A module over ii^T-_^(s/2) is then a diagonalizable f) -module W = ®^(^cW[fi], together 
with an L-operator L(z, A) G End(, {V W)) (a linear map commuting with h^^^ + h^"^^) 
depending meromorphically on 2;, A G C and obeying the relations 

(8) R'^^^\z-w, A-2r//i(3))L(i3)(^, A) L^^^\w, A - 2r//i«) = 

{z-w, A) . 
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For example, W = V, L{w,X) = R{w — zq, X) is a module over Er^r^ish), called the 
fundamental representation, with evaluation point zq. In more general examples of 
such modules were constructed: in particular, for any pair of complex numbers A, z 
we have an evaluation Verma module Mjy{z). It has a weight decomposition Ma = 
©°1o-^a[-A— 2j], with one-dimensional weight spaces Mj\,[fi]. The action of the L-operator 
is described explicitly in 0. Also, we have a notion of tensor products of modules 
over Er^nish)- The main examples considered in this paper will be tensor products 
AfAi(-Zi) ® ■ ■ ■ (S> MA„(-2n) of evaluation Verma modules and some of their subquotients. 

It will be convenient here to consider more general L-operators obeying the relations. 
So we define a functional module over £"^,^(5/2) to be given by a pair (W, L) where W is 
a space of complex-valued functions on a certain set and h acts on it as multiplication 
by a function, and A) is a meromorphic function of z and A acting as a difference 
operator on ® W , commuting with ® 1 + 1 ® /i, and obeying the relations (^. An 
example of such functional modules is provided by the "universal evaluation modules" 
of 1^. acts by multiplication by a continuous variable. Evaluation Verma modules are 
obtained by restricting the range of this continuous variable to a discrete set. 

For any module or functional module W over E-r^rjish), we define the associated op- 
erator algebra, an algebra of operators on the space Fun(iy) of meromorphic functions 
of A G C with values in W. It is generated by h, acting on the values, and operators 
a{z) , b{z) , c{z) , d{z) . Namely, let L{z) G End(V^ <S)Fnia{W)) be the operator defined by 
(L(z)(f (g) /))(A) = L{z, A)(f (g) /(A — 2rifi)) if hv = fiv. View L{z) as a 2 by 2 matrix 
with entries in End(Fun(Vr)): 

L{z){e[l]®f) = e[l]®a{z)f + e[-l]^c{z)f, 
L{z){e[-l](^f) = e[l]^b{z)f + e[-l](^d{z)f. 

The relations obeyed by these operators are described in detail in [[7| (in |^ these oper- 
ators are denoted by a{z), b{z) and so on). 

To each module we associate a central element of the operator algebra. It is given by 
the quantum determinant 

Det(z) = ^Ji^^(a{z + 2r^)d{z) - c{z + 2r/)6(z)). 

3.2. A class of representations by difference operators. Let zx^. . . , z„ e C be dis- 
tinct points and Ai, . . . , A„ G Z>o. Let us introduce difference operators acting on func- 
tions of n+1 complex variables A, xi, . . . , Xn- Let (T^^/)(A, xi, . . . , x„) = /(A, xi, . . . , Xi + 
a, . . . , x„) and (T^f){X, xi, . . . , x„) = /(A + a,xi, . . . , Xn). The steps a will always be 
±2t]. Let 

n n 

A_(z) = We{z - Zi + KiT]) and A+(z) = We{z - Zi - Kt]) . 

i=l i=l 
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The functions A±{z = —Xi), considered as multiplication operators, will be denoted 
simply by A±{—Xi). We also set s = J2i=ii^i + ^i)- 
With these conventions, we define: 



1=1 

n 



e{x) 

6{X + Z + Xi) -p-r 9{Z + Xj) 



- - E ^^^^ n ^A,(-., w. 



ciz 



^ e{-X + Z + X,-2s) yj e{z + Xj) ^ / ^ ^^.+2r,^.-2r, 



Det(z) =\Yd{z-Zi- AiT]) e{z - Zi + AiT] + 2r/) 



j=i 

Theorem 3.1. The difference operators a{z),h{z),c{z) together with d{z) defined im- 
plicitly by the determinant relation 

a{z + 2r])d{z) - c{z + 2T^)h{z) = ^%_M)Det(^) 

obey the relations of the elliptic quantum group (Qj with rjh = — Ylii^i + ^i)- 

Example. U n = 1, the generators act on functions of two variables Xi, A, and h acts 
as —Xi — Zi. If we introduce a new variable h = —r]~^{xi + Zi), so that the generator h 
acts by multiplication by h, we obtain a representation on functions of h, A. the action 
of the generators is given by the difference operators 

a{z)v{h,X) = 9{z - ^1 - r]h) ' f + ^'''^ v{h, A - 2r^) 

b{z)v{h,X) = + ' ~ ~ "^^^ 9{-r^h + A,r^)v{h + 2, A + 2r^), 

0{X) 

c{z)v{h,X) = JAzh±^_Z^l±^e{Tih + A,r^)v{h-2,X-2r^), 

d{z)v{h,X) = e{z-z^+ r]h) ~ "^f^ ~ ^'"^^ vih,X + 2r]). 

0{X) 

This is, up to normalization, the "universal evaluation module" of [0, Sect. 9. 



Proof of Theorem 15*. j| ; The proof consists of a straightforward verification of the sixteen 
relations. Let us give an example: one relation is 

, , e{z-w)e{X + 2r])^. . e{z - w - X)e{2r]) , , 
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This identity is verified by looking at each summand of b{z) and b{w) separately. The 
corresponding equation to one such summand typically looks like 

_ e{z - w)e{X + 2r]) d{X + w + Xk) .-i-r 6{w + Xi) 

e{z-w-27])e{x) e{\) ^f}j{xk~x,y 

X (11 + + - ''^'^ei!.^;t ^ ^ ^'^^ 

e{z-w- \)9{2r^) 9{X + ELil^. + ^^ + KyMw + Xk) , ^ 

+ e{z-w-2r,)eix) m (ll^(- + -^))x 

%\ei^^y eix-2rj) • 

By taking into account that each summand of the above equation involves a factor 

-j-j- e{z + Xi)e{w + g(A + Er=i(^» + + ^^^)) 

^(x.-x,) ^ ^(A) 

the task reduces to verifying 

e{X + w + Xk- 2ri) e{z - w)e{X + w + Xk)9{z + Xk- 2ri) 



9{z + Xk 



e{X - 2ri) e{z-w- 2r])9{X) 

e{z w X)0{2,i)0{a' + .r,)0{X + z + .r, 2r]) 



6{z-w -2r))6{X)6{X-2r)) 
which we write in the form 
fi{z, w, X) = f2{z, w, A) + fsiz, w, A) . 



This identity is proved in two steps. First, one shows that the functions fi{z,w,X), 
i — 1,2,3 , transform in the same way under A— >A + 1,A— >A + t,. The transformation 
laws thus obtained are the following: 

fi{z,w,X + l) = fi{z,w,X) 
Mz, w,X + r)^ e-2-(-+-'=)/,(z, w, A). 

for i = 1, 2, 3 . The second step is to show that the above relation holds for the residues 
of the functions fi{z,w,X). Here, one has to show that the identity holds for X — 2r] 
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and A = . For A = one obtains 








This proves 0, paragraph 3, relation 2. The other relations are proved in a similar but 
often more intricate fashion. 

One identity that is used is the vanishing of the sum of the residues at f = —Xj — 2//, 
—Xi of the function 



a consequence of the double periodicity f{v + l) = /{v + t) = f{v). 

Also some of the more complicated relations, such as d{z)d{w) = d{w)d{z) turn out 
to be consequences of the simpler relations and the fact that the determinant is central. 



Remark. This is the difference elliptic analogue of formulae that have appeared in the 
literature. In the rational difference and differential case such a formula has been written 



3.3. Restrictions. The linear difference operators defined in the previous subsections 
act on meromorphic funtions of complex variables A, xi, . . . , Xn- To compare these oper- 
ators with evaluation modules of the elliptic quantum groups and with transfer matrices 
of IRF models, we have to restrict their action to the space of meromorphic function 
defined on sub manifolds of C""*"^. The conditions for a difference operator X with mero- 
morphic coefficients to be defined on meromorphic functions on a submanifold 5* are 
that the value at a generic x E S of Xf{x) is well-defined (i.e., there are no poles at 
generic points of S) and is only a function of the values of / at points of S. Equivalenty, 
a difference operator X can be restricted to S if it maps functions vanishing on S to 
functions vanishing on S. The restriction is then identified with the induced action on 
the quotient by the function vanishing on S. 

These conditions are fulfilled in the following situations: 

1. Restriction to discrete values of Xi. We assume that Zi, rj are generic and that the 
Aj are non-negative integers. We take 5* to be the set 



6{2s + x^ + v + 2r]) yj 6{v - zi - Ai)9{v - + A; + 2r]) 
9{v + Xi + 2r]) Al 0(^y + xi)9{v + xi + 27]) 



□ 




5*0 = {{X,Xi,...,Xn)eC^'^\-Xi = Zi + r]{Ai-2mi), 
rrii = 0,1, ... ,Ai,i = 1, ... ,n}. 
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Since the steps in the difference operators are by muhiples of 2//, it is clear that one 
can restrict the action of a{z), . . . , d{z) to functions on subsets given by conditions 
Xi G ai + 2r]1,, for generic Oj, rj. The genericity condition on a^, r] is that the poles at 
Xi — Xj mod r of the coefficients of the difference operators are never on the subset. 
What we have to check is that the restriction to these finite sets of values for Xj is 
well-defined. Since only b{z), c{z) shift the value of Xi, it is sufficient to consider 
these two operators. For a function /(A,xi, . . . the value at —Xi = Zi + rjAi 
of b{z)f appears to depend on the value of / at —Xi + 2rj = Zi + rjAi + 2?7, which 
is not in 5*0, but in fact it does not, since the coefficient A^{—Xi) vanishes there. 
Similarly c{z) is well-defined on So- 

For this restriction h has discrete spectrum. It will be used to compare our 
representation with tensor products of irreducible representations. 
2. Restriction to X = rjh. Let S be the set 



5*1 = {(A, xi, . . . , x„) I A = - '^{xi + Zi)]. 



Then h{z),c{z) can be restricted to functions on Si. Indeed, if m is a function 
vanishing on 5*1, then T~^'^T'^u still vanishes on 5*1. The denominators of the 
coefficients of h{z) do not vanish at generic points of Si. Thus h{z)u = if m 
vanishes on Si. The same reasoning applies to c{z). 

This restriction is needed, as we shall see, to construct commuting transfer ma- 
trices. 

3. IRF restriction. Consider the restriction of 6(2;), c{z) on functions on S = SQf\ Si. 
If Zi and T] are generic, the only possible pole on S in the coefficients of these 
differential operators comes from the denominator ^(A). This denominator does 
not vanish if we assume for instance that the Aj are all odd. 

Here S is finite, so that the restriction is to a finite dimensional space of functions, 
which will be identified with the space of states of an IRF model. 

3.4. Commuting difference operators. One of the main properties of i?-matrices to 
statistical mechanics is that L-operators obeying quantum group relations give rise to 
commuting transfer matrices tryL. In |]lO| it is noticed that more generally one can 



consider try((i^' ^)L) for some endomorphism of V such that K ® K commutes with 
the i?-matrices. 

In our dynamical case it is known that the traces a{z) + d{z) commute for different 
values of z when acting on the zero weight space of a module over E^.^{sl2). Another 
possibility to obtain commuting operators is to take the twisted transfer matrix tiv{{K® 
1)L) with a suitable K: 

Proposition 3.2. For any & C, the operators T{z) = b{z) + i9c(z), 2; G C, restricted 
to functions on the suhmanifold Si given by the equation A + Yl^=i + Yl^=i = 0, 
form a commuting family. 

This proposition can be proved directly or by the following general argument. 
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It is first of all sufficient to consider the case = 1 since the other cases are obtained 
by conjugating the operator by the multiphcation by an exponential function of A. Then 
we may write T{z) as 



T{z) = J2 trvMiK(S)lL{z,X))T 
/i=±i 



with 

K 



1 

1 



The partial trace try[^] : End(V" C?> W) —>■ End{W) is the homomorphism such that if 
a G End(V), b G End(H^), then try[^](a ®b) = Y1 e\aei)h for any basis (cj) of V[ix\ with 
e* G V* defined by e*(ej) = 5ij and e*(w) = for w G V[u\, v ^ [i. We have 

(9) K®K R{z, A) = R{z, -A) K®K. 

We then write the RLL relations in the form 

z -w,X) ^ = 

L(2)(w,A)L«(2,A-2r//i(2)), 

multiply both sides by {K ® K)^^^'' from the left, and take the partial trace over a weight 
space {V ® Using (^), we obtain 

tr(v®y)M (^^''^ -w,-X + 2r,h)K^^^K^^^ 
L^^\z, \)L^^\w, A - 2r]hW)R(^^){z - w, X)-') = 
ti^v^vrnK^'^L^'Hw, X)K(^^L^'\z, X - 2r//i(2)). 

The next step is to use the cyclicity of the trace to bring the first R matrix to the right. 
For this we need the commutation relations of h with the product in the trace. Since 

+ h] = and h^'^K^'^ = -K^'^h^'\ we see that 

for any A G End(V^ V) commuting with h^-^^ + h^'^\ We then get 

R^'^^z - w, Xy'R^-'^\z -w,-X + 2r]{h + 2/i))T^^"'^ = 

T{w)T{z) 

If we then have a relation rihu{X) = Xu{X), and we apply the above equation to u, we 
may replace h in the left-hand side by X/t] — 2fi, and the R matrices cancel, so that 
T{z)T{w) = T{w)T{z). 
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3.5. Evaluation modules. Here we show that the restriction of the operators a, 6, c, d 
of Prop. |3.1| to functions on the submanifold 5*0 is essentially the tensor product of finite 
dimensional irreducible evaluation modules of 

Proposition 3.3. Suppose that rjjZi, . . . , Zn are generic. Let a, b, c, d be the difference 
operators defined in \3. ^ restricted to functions on Sq and let k,{z) = YYi=i (^{z—Zi—riKi)~^. 
Then a{z) = K,{z)a{z), . . . ,d{z) = K,{z)d{z) define an E^.f^{sl2) module isomorphic to the 
tensor product L\^{zi — r^) ■ ■ ■ (g) L\^{zn — rf) of irreducible evaluation modules. 

Let W be the space of functions on 5*. It is a vector space over the field of meromorphic 
functions of A of dimension nr=i(^« + To prove this proposition we have to identify 
a highest weight vector in t>, i.e., an eigenvector of a{z), d{z) and h killed by c{z). 
The eigenvalues {A{z, A), D{z, A), A) of {a{z), d{z), h) determine then by ^ uniquely an 
irreducible module up to isomorphism. 

Let Sa{xi) E W he the delta function at Xi = a: Sa{xi) = 1 if Xj = a, Sa{xi) = 0, if 
Xi 7^ a. The highest weight vector may be taken as the product of delta functions 



Y[6^,^_A^r,{Xi) . 



i=l 



This function is indeed annhilated by c{z) and hvh.w. = JZ-^i'^h.w.- Moreover Vh.w. is an 
eigenvector for a{z) and d{z): 



= n A.,) ^(A^Sl|^^±SlM,„,. 

n 

YlOiz - Zi- K'q)vh.w. ■ 



1=1 



Thus a{z)vh.w. = A{z, X)vh.w., with A{z, A) = nr=i ^(^ - - ^iV)- Similary d{z)vh.w. 
D{z, X)vh.w. with eigenvalue 

^(^'^) = ~ !!!F^' ^'^ Det{z - 277, - 2r/, A + 2r/) . 

^(A - 2r] Z7=i ^i) A ^iV - 2r7)g(^ - Zj + AjT]) 

e{X) Al e{z-z,-A,r]-2r]) 

Thus the eigenvalues of {a{z), d{z), h) are (1, D{z, A), ^ Aj) with 

,,_ ^(A-2r/Er=iA.) fr g(^-^. + A.77) 
m fj e{z -z.- A.^) ' 



which indeed reproduces the highest weight defined in p. 750. 
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3.6. Separation of variables: continuous case. In this section we find an analogue 
of the results of ^]^, |2.5| . We consider the continuous case, in which the variables Xi take 
arbitrary complex values. Then the eigenvalue problem for a function u on 5*1 reads 
T{z)u{x) = e{z)u{x), with 

6{\ — z + Xi) Yj 9{z — Xj) 



T{-z)u{Xi,...,Xn) = WV\ 



,=1 W j^,0{x.-x,] 



n 



( JJ Oixi + Zk + r]Ak)u{xi, ...,Xi-2T],--- , x„) 

k=l 
n 

+ Y\_ ^i^i + Zk- r]Ak)u{xi, ...,Xi + 2r],..., , 

k=l 

where we view a function u on Si as a function of Xi, . . . , x„ by setting A = — "^{xj + Zj). 

From this formula it is clear that T{z)u{x) is an entire holomorphic function of z with 
theta function behavior as z is translated by elements of the lattice F = Z + rZ. It 
follows that a necessary condition of e{z) to be an eigenvalue is that e belong to the 
space Thetanixo) of theta functions of order n with character xo : T — C* (see the 
Appendix) given by Xo(^ + st) = (— 1)"'^^"'"^) exp(27ri ^ 2;^). This means that e is an 
entire function obeying 

e{z + 1) = xoimz), e{z + r) = xo(r)e-™(2^+-)e(^). 

The method of separation of variables consists in looking for common eigenfunctions u{x) 
in the factorized form u{x) = Yli=iQ{xi). Setting z = —Xi in the eigenvalue problem 
(T(z) — e{z)) u = we see that a necessary condition is that Q,e obey the difference 
equation 

n 

A+{x)Q{x-2t]) + A_{x)Q{x + 2r]) = e{-x)Q{x), A±{x) = JJ ^(x + ± r^A^) 

k=l 

As explained in the Appendix, this difference equation has an elliptic polynomial solu- 
tions, i.e., a solution of the form 

m 

(10) Q{x)=e^''l[9{x-Wk). 

k=l 

if ^ Aj is an even integer 2m. Such a solution may be constructed by the Bethe ansatz: 

Proposition 3.4. Suppose that Ai + ■ ■ ■ + A„ = 2m for some positive integer m, and 
let (a, wi, . . . , Wm) he a solution of the system of Bethe ansatz equations 

(11) A+{wi) JJ e{wi - Wj - 2r]) = e'^^'^A.iwi) JJ 9{wi -wj + 2r]) i = l,...,m, 

such that Wi 7^ Wj mod F, (i 7^ j). Then ^ = 11 Qi^i) ''^'^^h Q{x) = e""^ YYk=i ^(^ ~ "^fc) 
is a common eigenfunction ofT{z). 
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Proof: This is a rephrasing of the first part of Prop. [A .41 . □ 



Definition. An eigenfunction of the form of Prop. |3.4| is called Bethe eigenfunction. 

Conversely, let us suppose that ^ Aj = 2m, m G Z>o and show that all eigenfunctions 
in a suitable class are of this form. Let, for a character x : F — »• C*, Hmix) be the space 
of meromorphic functions of n complex variables Xi, . . . , x„ such that 

«(• • • ,Xi + 1, ■ ■ ■) = ■ ■ ■ ■ ■), 

The following result can then easily be verified using the behavior of the coefficients of 
the difference operator T{z). 

Lemma 3.5. For any character x o-nd any 2; G C, T{z) preserves 'Hm{x)- 

Let for a character x G F*, C Om(xo) be the set of functions e so that there 

exists a holomorphic common eigenfunction u G Hmix) of T{z), z G C, with eigenvalue 
e{z). 

Theorem 3.6. Suppose that Ai + ■ ■ ■ + A„ = 2m for some positive integer m. Then 
e G if and only if there is an eigenfunction u{x) of the form u{x) = YYi=iQ{^i) 

with eigenvalue e, such that Q{x) = e"-^ YYiLi ^(^ ~ '^i) /^'^ some solution (a, wi, . . . , Wm) 
of the Bethe ansatz equations with 

X{1) = (-l)™e^ xir) = (_i)'»e"-+2-iE«'fc. 

Proof: It remains to show that if e G then this eigenvalue corresponds to a Bethe 

eigenfunction. 

Suppose that e G and m is a holomorphic eigenfunction in Timix) "with this 

eigenvalue. Then, for each i, the function u viewed as a function of Xj belongs to Q>m{x)- 
By setting z = —Xi in the eigenvalue equation T{z)u{x) = u{x), we see that m is a 
solution of the separated equation 

n 

(12) Yl e{xi + Zk + v) T-^^'uix) + e{xi + Zk-7i) Tly{x) = e{-Xi)u{x). 

k=l 

Thus this equation admits a non-trivial solution Q{xi) in Qmix) (we consider the re- 
maining variables Xj, j ^ i as fixed). By the factorization theorem for theta func- 
tion (see the Appendix), such a solution is, up to normalization, of the form Q{xi) = 
e"^ (^{xi — Wj), and x is related to a and Wk by the equations stated in the Theorem. 
Setting Xi = Wj in (|12|) gives then the Bethe ansatz equations. □ 
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4. IRF MODELS WITH ANTIPERIODIC BOUNDARY CONDITIONS 

We consider here the case where Ai = ■ • ■ = A„ = 1 and show that the commut- 
ing transfer matrices T{z) restricted to functions on Sq H Si are transfer matrices of 
IRF models with antiperiodic boundary conditions. The IRF (interaction-round-a-face) 
models (see 0, |T|, |^) of statistical mechanics are two-dimensional lattice models in 
which the configurations over which one sums in the partition function associate an 
element, the height, of a certain set to each pair of neighboring points. The weight of 
a configuration is the product of local Boltzmann weights associated to each "face" , or 
square formed by four neighboring points. The local Boltzmann weight W{a,b,c,d\z) 
depends on the heights a, b, c, d at the neighboring points and the spectral parameter 
z e C. The spectral parameter can be any fixed number, but more generally, in the 
inhomogeneous model, one associates a spectral parameter to each row and column of 
faces in the lattice and takes z to be the difference between the row parameter and the 
column parameter. In the simplest s/2 case the heights are real numbers and the local 
Boltzmann weights are related to matrix elements of our dynamical i?-matrix by 

(13) R{z, X = -2r]d) e[c - rf] (g) e[6 - c] = ^ W{c, b, a, d\z) e[b - a] (g) e[a - d], 

a 

where c — d,b — c,b — a,a ~ d & {^l? If latter conditions are not fulfilled, then 

w = o.g 

As shown in [^, the transfer matrix a{z) + d{z) on the zero weight subspace of the 
tensor product of two dimensional representations is then identified with the row-to-row 
transfer matrix of the IRF model with periodic boundary conditions. Here we show 
that a similar computation applied to the twisted transfer matrix b{z) + c{z) gives the 
row-to-row transfer matrix of the IRF model with antiperiodic boundary conditions. 

The formulae are as follows: Let V{zi) ® ■ (g V{zn) be a tensor product of two di- 
mensional evaluation modules and assume that n is odd. Then the L-operator for this 
module is 

n 

L{z, A) = R^^^\z -zi,X-2r]J2 h^^) ' • -^^"""^H^ - ^n-i, A - 2r]h^''^)R^°''\z - A). 

1=2 

Then the twisted transfer matrices, see acting on functions m(A) restricted to the 
submanifold given by the equation X = rjh form a commuting family parametrized by the 
spectral parameter. Let us call these transfer matrices Tirf(-z). The restriction means 
that Tirf(-2) preserves functions u{X) of the form u{X) = J2^_l^r|^^W'^[^A some vectors 
u[fi] of weight /i. Here Sa{b) is one if a = 6 and zero otherwise. The restriction to this 
subspace of function of A is well-defined if n is odd, as the singularities of the L operator 
are at A G 2riZ whereas A takes values rjii, where /x, the weight of a vector in ^V{zi) is 
an odd number if n is odd. So if rj is generic, there are no singularities at this points. 



^The normalization chosen here is so that W{1, l + l,l + 2,l + l\z) = 1 rather than the more common 
6{z — 2r])/9{2ri). Also our spectral parameter is normalized in a different way than in the literature. 



20 



Explicitly, 

Tirf{z) = h{z) + c{z) = 

(14) - A - 2?7/i('^))i?(0")(^ - A)r-'^'^, 

with 

X = 



1 

1 

The product of i?-matrices acts onV ® {^'^^-^V), where the index (0) refers to action on 
the first factor of the tensor product, tr^^^ denotes the trace over the first factor. The 
subsequent factors are numbered accordingly from 1 to n. 

The twisted transfer matrices are defined on a 2"-dimesnional space. This space has a 
basis J] (T,(^)e[o"i] eg) ■ ■ ■ e[o'„] with (jj e {1, — 1}. It is convenient to write this basis in 
terms of antiperiodic paths: Let for ai, . . . , a„, a„+i G Z + n/2, such that \ai — a^+il = 1 
(i = 1, . . . , n), and a„+i = — ai 

|ai, . . . , an+i) = (^-27,o„+i(A)e[ai - 02] (8) • • • (8) e[a„ - a„+i] 

Proposition 4.1. For any antiperiodic path |ai, . . . , a„, a„+i = — oi); w^e /lave 
Tikf{z) |ai, . . . , a„+i) = 

n 

(15) ^ JJiy(aj+i,ai,6j,6i+i) . . . , , 

,6„,6„+i=— 61 i=l 



This expression is, by definition, the transfer matrix of IRF models with antiperi- 
odic boundary conditions. The partition function for the IRF model with antiperiodic 
boundary conditions in one direction is then trT(-u;i) ■ ■ ■T{wm)- the Wi are the spectral 
parameters associated to the m rows, and the Zi are the spectral parameters associated 
to the n columns. 

Proof of the proposition: Let e*[l], e*[— 1] be the dual basis to e[l], e[— 1]. 

Tirf{z) |oH, ■ ■ ■ , On+l) — 

n n 

J2 t4l]K^°^R'-''^ {z-z,,X-2r)J2 h^'^)R^'''^ {z - ^2, A - 277 J] /i») 

/u i=2 i=3 

■ ■ ■ R^'^''\z - Zn, A) 5_2r,(a„+i-/i)(A)e[ai - 02] (8) • • • (8) e[an - On+l] = 
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e^°^[/x] (g) e[ai - 02] ® • • • ® e[a„ - a„+i] = 

n 

(^-2T;6„+i(A)e(°)[a„+i - hn+i] ® e[ai - 02] (g) • • • e[a„ - a„+i] = 

i^(On-i)(^ - A - 277/i(-)) W^(a„+i, a„, &„+i)(5_2,6„+,(A) 

e^°^ [an - hn\ <S> e[ai - 02] <8) • • • ® e[an-i - a„] e[6„ - 6^+1] = • • • = 

n 

J2 *K+i - - ^1, A - 2r7 5^ • • • 

6„+i,...6i 1=2 

n 

j=i 

5-2„&„+i(A)e(°)[ai-6i]® 
e[ai - 02] ® e[aj_i - Oj] (g) e[6j - • • • e[6„ - = • • • = 

6„+i,... ,61 

5-2^fe„+i(A)e(°)[ai - hi] ® e[hi - 62] ® ■ ■ • ® e[6„ - 6^+1] = 
^ e(°)*[a„+i - 6„+i]e(°)[-ai + 61] ® e[W - 62] <8) • • • e[6n - K+i] 

6„+i,... ,61 

(A) n«^,W(a,+i 



By evaluating the linear form e*[a„+i — 6„+i] we see that only the terms with hn+i = —hi 
contribute to the sum. The proof is complete. 



4.1. Separation of variables for IRF models. Let us consider the eigenvalue prob- 
lem for the transfer matrix restricted to 5*0 fl 5*1 for A 1 = • • • = A„ = 1 and n 
odd. The eigenfunction may be viewed as a function u{xi, . . . ,Xn) defined for Xi e 
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{—Zi — Tj, —Zi + 7]}, 1 <i <n. The eigenvalue problem reads T{z)u{x) = e{z)u{x), with 

>A e(A - Z + Xi) -p-r ^(2; - Xj) 

T(-.Mx„...,x„) = \}^0i^^) 

n 

(16) ( JJ 9{xi + Zk + r]Ak)u{xi, ■ ■ ■ , Xt - 2r], ■ ■ ■ , Xn) 

k=l 

n 

+ Y\_ ^i^i + Zk- r]^k)u{xi, ...,Xi + 2r],..., Xn)) , 
fc=i 

where A = - XlLi(^fc + Zk)- 

Theorem 4.2. Suppose that Ai = ■ • • = A„ = 1 with n odd, 1] ^T, and that Zi 7^ Zj+2r]£ 
mod r for i ^ j and £ = 0, ±1. Let T{z) = h{z) + c{z) he the transfer matrix restricted 
to the 2" dimensional space of functions on So fl ^i. Then a function e{z) is a common 
eigenvalue of the transfer matrices T{z), z & C, if and only if 

(i) e e Bnix) with = (-1)", = (-l)V^^S^^ and 

(ii) e obeys the quadratic relations 

n 

(17) e{zi-r])e{zi + r]) = Y\_0{zk- Zi + 2r])9{zk- Zi-2r]), i = l,...,n. 

k=l 

To prove this theorem, let 5* = x^=i{—Zi — rj, —zi + r/}. Let us first assume that u{x), 
X G 5 is a common eigenfunction of T{z) with eigenvalue e{z). In particular u does not 
vanish identically on S. From the transformation properties of T{z) under shifts of z 
by r we see that e{z) has to belong to 0n(x)- Then setting z = —Xi in the eigenvector 
equation T{z)u{x) = e{z)u{x), we get the separated equation 

n 

(18) Yl 9{xi + Zk + v) T-^M^) + G{xi + Zk-v) Tly{x) = e{-Xi)u{x). 

k=l 

Inserting the two values of Xi yields 

(19) ^ 

riLi ^i^k - Zi + 2r]) u{xi, . . . ,-Zi - rj, . . . ,Xn) = e{zi - r])u{xi, . . . , -Zi + r], . . . ,Xn), 

riLi ^(^fc - Zi- 2r]) u{xi, ...,-Zi + r],...,Xn) = t{zi + ?7)m(xi, . . . , -2;^ - 77, . . . , x„). 

By multiplying these two equations we see that e must obey the identity ([TtI), provided 
we can prove that, for at least one choice of Xj, the product 

m(xi, . . . , -Zj - ?7, . . . , x„) m(xi, . . . , + ?7, . . . , x„) 

does not vanish. This follows from the fact that u is not identically zero, so that at least 
one factor of this product is nonzero, and that the product of theta functions on the 
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left-hand side of (|T9D is not zero with our assumption on the zj, so that also the other 
factor is nonzero. 

We have thus shown that a necessary condition for a funtion e{z) to be a common 
eigenvalue is that e is a theta funtion obeying the quadratic relations (p!7|). 

Conversely, let us suppose that e obeys (|T^. Then, for every i, the system of equations 

n 

Y\ d{x + Zk + r]) Q{x - 27]) + 9{x + z^ - r]) Q{x + 2r]) = e{-x)Q{x), x = -Zi± t], 

k=l 

admits a non-trivial solution Qi{x). It follows that u{x) = YYi=i Qii^i) obeys (|TB|). Thus, 
for any x ^ S, {T{z) — e{z))u{x), viewed as a function of z is a theta function in 0„(x) 
vanishing at n distinct points —Xi, . . . , —Xn- As explained in the Appendix, this implies 
that that either the function vanishes identically or — Xl^j = Xl-^*- Since for x ^ S, 
^ Xj is —Y^Zi plus an odd multiple of t], the latter alternative cannot hold. Thus u is 
an eigenfuntion with eigenvalue e. The Theorem is proven. 



5. Conclusions 

In this paper we have studied integrable models associated to elliptic curves by 
Sklyanin's method of separation of variables. We have shown how in the elliptic version 
of the Gaudin model, previously considered in 0] , this method implies the completeness 
of Bethe eigenfunction in Verma modules, in the sense that for every eigenvalue one 
has a Bethe eigenfunction (although in the case of degenerate eigenvalues one cannot 
exclude the existence of additional eigenfunction not of Bethe type). 

In the difference case, the analogue of the transfer matrix in the separated variables 
was shown to be a twisted transfer matrix associated to a representation of £',-,77(^/2) by 
difference operators. The eigenvalue problem can be posed in two different cases. In the 
continuous case, the variables are assumed to be complex and one asks for eigenfunctions 
with theta functions properties, as one does in the differential case. These eigenfunctions 
are of Bethe ansatz type. In the discrete case, the variables are assumed to take values in 
a finite set. The eigenvalue problem can still be solved by separating variables yielding 
eigenfunctions of the (inhomogeneous) row-to-row transfer matrix of IRF models with 
antiperiodic boundary conditions. In both cases one has a completeness result. 

Let us conclude by mentioning some open questions. In the difference case one has the 
transfer matrix in the "separated variables" . Is there a "quantum Radon transform" as 
in the differential case, which maps the transfer matrix of the IRF model with periodic 
boundary conditions to this transfer matrix? Is there an analogue of the equivalence 
between the local problem and the global problem of |T2[ in the elliptic case? 



Appendix A. Elliptic polynomials 

A.l. Theta functions. Let Imr > and set q = e^'"*'^. Let T be the lattice Z + rZ and 
r* ~ (C^)^ the group of group homomorphisms F C^. Let be the homomorphism 
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^'■X^ 2^(ln x(t) -rln x(l)) (or, more invariantly, 2^(a;iln xM), for 

any oriented basis uji,uj2 of F) from V* to the elliptic curve E = C/F. 

For X G F* let Qk{x) be the space of theta functions of level k and character x- 
It consists of entire holomorphic functions f{z) such that f{z + r + sr) = x(r + 
sr) exp(— 7riA;(s^T + 2sz))f{z) for all r + sr e F. 

The dimension of Qk{x) is zero ii k < 0. It is A; if /c > 1. The dimension of ©o(x) is 
one if 0(x) = and vanishes otherwise. 

We have the unique factorization result: 

Proposition A.l. The function of z e C 

k 

f{a,w;z) = e^^l[e{z-wj) 
j=i 

belongs to Qkix), with xir + sr) = (^_iYr+^)k^ra+s{aT+27:iEjWj) _ Every function in Qkix) 
is of the form C ■ f{a,w; z) for some constant C and this representation is unique up 
to permutation of the Wj if one requires the Wj to be in the fundamental domain F — 
{x + yT\x,y e [0,1)}. 

Proof: It follows from the transformation properties of theta functions that the number 
of zeros {2ni)~^ jQp,dlng, counted with multiplicities, in F of a theta function g G 0fc(x) 
is k. li Wi, . . . ,Wk denote the zeros of g then g{z)/ f{a,w; z) is doubly periodic and 
regular, thus constant. Uniqueness follows from the fact that a is uniquely determined 
by the zeros Wj and the character x- CH 

Corollary A. 2. Let E be the elliptic curve C/F, and, for k = 1,2, ... , S''{E) = E/Sk 
he its k^^ symmetric power. The map P(Ofc(x)) S^{E) sending a function to the set 
of its zeros modulo T, is injective. Its image consists of classes [wi, . . . ,Wk] such that 
X^j- Wj — (f){x) + kS, where S is the image in E of {1 + t)/2. 

A. 2. Interpolation. 

Theorem A. 3. Suppose zi, . . . ,Zk G C are pairwise distinct modulo F and x £ F* is 
such that Y^i=i 7^ 4>{x)'^k5 mod F. Then for any yi, . . . , G C there exists a unique 
function f G ©fe(x) ^^c/i that f{zi) — yi for all i — 1, . . . ,k. 

The interpolation formula giving / is 
Here a and b are such that the character of all terms in the sum are equal to x, so 



K 



1+T 



25 



for any choice of the branch of the logarithm. The assumption on x ensures that h 
so that the denominator 9{h) does not vanish. 

The function is unique since the difference of any two is a theta function vanishing at 
m points. By Corollary [A.2| , with our assumption on x, it must vanish identically. 

A. 3. Difference equations. Consider the following problem arising in integrable mod- 
els. Given 7 G C — F and functions A^{z) G 0fc(x+), ^-{.z) G Qk{x~)i S^id e{z) such 
that the difference equation 

(20) A^{z)Q{z - 7) + A_{z)Q{z + 7) = e{z)Q{z), 

has a non-trivial solution Q{z) in some 9m(x)- 

A necessary condition is that all terms are theta functions with the same character. 
So e has to be of level k and = X-i^)-, X+i^')^^'^^"'"^ = X-{'^)^~'^'^^'^"^ = character 

of e. 

A pair (e, Q) of theta functions (with Q non-trivial) obeying the difference equation 
(^) will be called an elliptic polynomial solution of (PPI). 

Proposition A. 4. Suppose A±{z) G Ofc(Xd=) with x+('" + st) = X-(r -|- sr)e~^'^*''"™'^ for 
all r + ST & r and some positive integer m. If a,wi, . . . , Wm are solutions of the system 
of equations 



(21) A+(m,) Yl 9(10, - w, - 7) = e^'^^A^w,) J] 0(m, - w, + 7) 
such that Wi 7^ Wj mod F for all i j , then the functions 



m. 



(22) Q(z) = e-n^O 



and 

A+iz)Qiz-^) + A4z)Qiz + ^) 



(23) e{z) 



Q{z) 



form an elliptic polynomial solution of (\2q). Conversely, if{e, Q) is an elliptic polynomial 
solution of ^dj), then there exists a solution a,Wi, . . . ,Wm of the system such that 
Q is of the form ^2^ ) (up to a multiplicative constant) and e is given by (^^. 

Proof: Let Q be the function defined in (0). The ratio ( p3D obeys the transformation 
property of a theta function, but may be singular at the zeros Wi of Q. 

The i^^ equation in the system, or more precisely, the equation equivalent to it if 
7 G C-F: 

m m 
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is the condition that the left-hand side of the (pOD vanishes at Wi. Thus, if a, ti^i, . . . , Wm 
obey the system of equations, then the quotient e of this left-hand side by Q{z) is regular 
at Wi and thus everywhere and (e, Q) is an elliptic polynomial solution. 

Suppose now that (e, Q) is an elliptic polynomial solution. By Proposition |A.1| , Q 



can be written in the form ( p2D if we normalize it properly. The points Wi are the zeros 
(modulo r) of Q. Then the left-hand side of the difference equation vanishes at Wi so 
the Wi are a solution of (|2T| ) □ 

Remark. The system of equations is usually written in the form 

ejWi - Wj - 7) _ Oa^A_iWi)_ 
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